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Abstract 

The restricted three-body problem describes the motion of a massless particle under the 
influence of two primaries of masses 1 — \i and /x that circle each other with period equal to 
2tt. For small fi, a resonant periodic motion of the massless particle in the rotating frame can 
be described by relatively prime integers p and q, if its period around the heavier primary is 
approximately 2-Kp/ q, and by its approximate eccentricity e. We give a method for the formal 
development of the stable and unstable manifolds associated with these resonant motions. We 
prove the validity of this formal development and the existence of homoclinic points in the 
resonant region. 

In the study of the Kirkwood gaps in the asteroid belt, the separatrices of the averaged 
equations of the restricted three-body problem are commonly used to derive analytical approx- 
imations to the boundaries of the resonances. We use the unaveraged equations to find values 
of asteroid eccentricity below which these approximations will not hold for the Kirkwood gaps 
with q/p equal to 2/1, 7/3, 5/2, 3/1, and 4/1. 

Another application is to the existence of asymmetric librations in the exterior resonances. 
We give values of asteroid eccentricity below which asymmetric librations will not exist for the 
1/7, 1/6, 1/5, 1/4, 1/3, and 1/2 resonances for any \i however small. But if the eccentricity 
exceeds these thresholds, asymmetric librations will exist for fi small enough in the unaveraged 
restricted three-body problem. 

1 Introduction 

The restricted three-body problem describes the motion of a massless particle under the influence 
of two primaries of masses 1 — fi and /j. The Hamiltonian 

H = ¥ pl + * + yp * - xp y ~ (J+yiy/i ~ » (((x- 1)2 + ^1/2 - *) (i-i) 

gives the equations of motion of the massless particle. In (|1.1|) . it is assumed that the primary of 
mass 1 — /x is located at (0,0), that the primary of mass \i is located at (1,0), and that the frame 
of reference rotates with the second primary in the anticlockwise sense with period 2ir. 
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In terms of the heliocentric Delaunay variables L,l,G,g, the Hamiltonian becomes 



H = -G- »n'(L, G ,l,g) = ~-G- tfl{L, G, I, g) - (1.2) 

where 

O' = —7- — r cos 9 

(1 + r 2 - 2r costf) 1 / 2 

n = n'-i/L 2 . (1.3) 

The variables r, 9 can be obtained in terms of L,l,G,g using the equations 

e = (1 - G 2 /L 2 ) 1/2 

a = L 2 /(l-^) 

I = E — e sin £7 
cos f = (cos E — e)/(l — e cos 
sin 1/ = (1 - e 2 ) 1/2 sin E/(l - e cos E) 

r = a(l - ecosE). (1.4) 

The Cartesian coordinates x,y in the rotating frame used in (jl.lj) are given by x = rcos9 and 
y = rs'm9. When [i = 0, the orbit of the massless particle in the inertial frame is an ellipse with 
eccentricity e and semimajor axis a; I, v, and E are the mean, true, and eccentric anomalies; g 
is the argument of the perihelion in the rotating frame and G denotes angular momentum. The 
variables r and 9 are the polar coordinates of the massless particle in the rotating frame. 

We use L, I, G, g to investigate motion near resonances as the Hamilton's equations of (|1.2|) 
take an especially simple form for /i small. If L = (p/q) 1 ^ 3 , G = (p/o , ) 1 ^ 3 (l — e 2 ) 1 / 2 , < e < 1, 
and [i = 0, the massless particle moves on an ellipse with period 2irp/q and eccentricity e in the 
inertial frame. It is assumed that p and q are relatively prime positive integers. The motion is 
periodic in the rotating frame as well with period equal to 2irp. If I and g are multiples of ir and the 
unperturbed orbit does not collide with the second primary, these periodic motions in the rotating 
frame perturb to periodic motions for /i > and small pQ. Given p, q, e it might appear that there 
are four possibilities as I and g can be either or ir, but in fact only two of these are distinct. These 
are the two q/p resonant periodic motions of the restricted three-body problem that correspond to 
eccentricity e. 

The Hamiltonian H is conserved by the flow, and for a given H and fi small it is possible to 
solve (|1.2|) for G using the implicit function theorem. Thus the Poincare section g = can be 
identified with the l-L plane. Each resonant periodic motion corresponding to p, q, e intersects this 
Poincare section exactly p times. Typically, if one of the resonant periodic motions is of elliptic 
type the other is of hyperbolic type jllj . In Figure ^ we have shown the periodic points on the L-l 
section for q/p = 3/1 and q/p = 7/2 with elliptic points marked as circles and hyperbolic points as 
crosses. The perturbing term Q and the 0(fj, 2 ) term in (|l,2j) are unchanged by the transformation 

L <— L, I < 1, G <— G, g < g, which is symplectic with multiplier —1. This discrete symmetry 

of the Hamiltonian (|1.2j) has the following implication for the return map to the l-L section: if 
(Iq,Lq) maps to (h,Li) then the return map sends (— h,L\) to (— Iq,Lq). 
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Figure 1: The l-L plane shown above is the Poincare section g = 0. Periodic points as well as stable 
and unstable manifolds are shown for q/p = 7/2 and q/p = 3/1. 

We obtain a scaled version of the return map near q/p resonance in Sections 2 and 3. The 
stable and unstable manifolds of the hyperbolic points of such a return map nearly coincide and 
the angle of transversality can be upper bounded by a quantity that is exponentially small in the 
small parameter fx. A discussion of exponential splitting of separatrices can be found in the work of 
Gelfreich and Lazutkin [3], Fontich and Simo [I], and Holmes, Marsden and Scheurle [7j. Instability 
is often associated with resonance, and it is therefore natural to look for transverse homoclinic points 
near resonances. However, it appears that there has been no construction of transverse homoclinic 
points near resonances in either the restricted three-body problem or in some other version of the 
planetary problem. The rescaling of the return map given in Section 3 brings this problem into 
sharper focus and it is possible to make an analogy to the discussion of the standard map in p]. 

In Section 4, we describe a procedure for the formal development of the stable manifold of a 
resonant periodic solution in powers of fi 1 / 2 . This procedure is specially adapted to the restricted 
three-body problem. In Section 5, we have included a proof of validity of this formal expansion. 
We also prove the existence of homoclinic points. A part of the verification essential for the proof 
is carried out in Section 6. 

Although investigations of the Kirkwood gaps in the asteroid belt have used physical models 
that include the secular variation of Jupiter's elements and the effect of Saturn, the averaged 
circular restricted three-body problem is still used to approximate the boundaries of the resonant 
regions. The more complicated models are essential to explain the dynamics within the resonant 
regions. An account of these models and their use can be found in the monograph by Morbidelli 
jH]. The use of the averaged equations of the circular restricted problem to sketch the boundaries of 
resonance can be found in Dermott and Murray [Sj , Henrard and Lemaitre |S] , and Lemaitre jS] . It 
appears to be known that the boundaries obtained from the averaged circular restricted problem do 
not work well at low eccentricities [H] ^3]. In fact for the q/p = 2/1 case, the left boundary cannot 
even be computed near e = jHj. In Section 7, we interpret the calculation of the boundaries in 
terms of the unaveraged circular restricted problem. For the commonly studied Kirkwood gaps, we 
give values of the eccentricity below which the approximation of the boundary will not be valid. 

Asymmetric librations near exterior resonances with p > q in the averaged circular restricted 
problem have been investigated by Beauge [2| and more recently by Voyatis et al. ^2]- I n Section 
8, we show the existence of these librations in the unaveraged circular restricted problem for q/p 
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equal to 1/2, 1/3, 1/4, 1/5, 1/6, and 1/7. We also give minimum values of eccentricity for each of 
these exterior resonances which must be exceeded for the asymmetric librations to exist. 



2 Pert ur bat ive form of the return map 

Let L range over the interval [(p/q) 1 ^ 3 — Sl, (p/q) 1 ^ 3 + $l] for some 5l > 0. We assume p/q ^ 1/1. 
Let 6 range over the interval ^rnin^^-max 

] such that L 2 (l — e) > 1 + 5 if p > q and such that 
L 2 (l + e) < 1 — 5 if p < q, for all allowed values of L and e and some 5 > 0. In addition, assume 
< e m j n and e max < 1. We take the range of the Hamiltonian H to be the set of values of 
— (p/q)~ 2 ^ 3 /2 — (p/q) 1 ^ 3 ^ — e 2 ) for e m i n < e < e max . Then for any allowed value of H and L, 
any real values of the angles I and g measured modulo 2tt, and \i sufficiently small, the implicit 
function theorem enables us to solve (jl.2|) for G uniquely. In fact, G will be an analytic function 
of l,L,g,H,ix. 

We use either g = or g = it to define the Poincare section. Since Hamilton's equations formed 
using the Hamiltonian (|1.2j) imply that dg/dt = — 1 + O(fi), the return map is well defined for L 
and H in the intervals specified by the previous paragraph, any value of I, and // small. Since H 
is conserved by the flow, we may identify the Poincare section for fixed H with the l-L plane as in 
Figure^ The return map preserves the area element dl dL. The return map will be denoted by T\ 
and Tf will be denoted by T p . In this section, we will obtain the perturbative form of T p . 

The Hamilton's equations of (|1.2|) imply that I = 1/L 3 — + 0(/i 2 ), g = —1 — /ufic, and L = 
1x0,1. We seek a solution of these equations with the initial conditions 1(0) = Iq, L(0) = Lq, g(0) = go, 
where go is either or it, and G(0) = Gq. It is understood that G must be obtained by solving 
(|1.2jl . If the solution is represented as l(t) = l a (t) + h(t)n + 0(fi 2 ), g(t) = g a (t) + <#>(i)/i + 0(n 2 ), 
and L{t) = L a (t) + L b (t)fi + 0(/u 2 ), then l a (t) = l + tjL%, g a (t) = g - t, and L a (t) = L . In 
addition, l^, and Lf, must satisfy 

h = (—3/La)L b — Ql, 9b = -^G, L b = Q t , 

where the partial derivatives of Q must be evaluated at I = Iq + t/L^, g = go — t, L = Lq, and 
G = Gq. By solving the equations above, we get 

l(t)=l + t/Ll + n\- f SlLdt-Tif I FtidtdT) +0{fi 2 ) 
\ Jo L o Jo Jo J 

g(t) = go-t-fi [ n G dt + o{fi 2 ) 

Jo 



L(t)=Lo + n nidt + 0{fi% (2.1) 
Jo 

where the partial derivatives of Q must be evaluated at I = Iq + t/L^, g = go — t, L = Lo, and 
G = Gq. The solution given by (|2.1j) is valid over any finite interval of time for initial conditions 
in the domain already indicated and for [i sufficiently small. 

To approximate T p , it is necessary to find the time t r at which g(t) = —2-Trp + go. Using the 
equation for g(t) in (|2.1|) . we get t r = 2-Kp — \x J Q 27rp dt + 0(fi 2 ). Using ()2.1|) . we may deduce that 
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Tp is given by 



h = h + 2irp/L 3 + m( —73 / n G dt- n L dt--j 



2irp 

2\ 



L\ = L + fi / Qidt + 0(n 
Jo 




where the partial derivatives of Q must be evaluated at I = Iq + t/L^, g = g$ — t, L = Lq, and 
G = Gq. The expression for T p given by IJ2.1J) is valid for L G [(p/^) 1//3 — <5l ; (p/q) 1 / 3 + <^i], for any 
real value of the angle Z, for H within a range that ensures avoidance of collision with the second 
primary as specified earlier, and for fj, sufficiently small. Since this domain of validity is compact, 
the 0(fj?) terms in 1)2. 2 Jl hold uniformly over the domain. 

The lemmas below are related to the first return map T\ and its pth iterate T p . 

Lemma 2.1. Assume that (Iq,Lq) maps to (l\,L\) under T\ (or T p ). Then (— li,Li) maps to 
(—lo,Lo) under Ti (orT p ). 

Proof. The Hamiltonian (|1.2|) is unchanged by the transformation / < 1, L <— L, g < g, G <— G. 

Therefore, if l(t), g(t), L(t), G(t) is a solution of the Hamilton's equations of Q1.2[) for < t < t* , 
then —l(—t), —g(—t), L(—t), G(—t), where — t* < t < 0, is also a solution. The lemma follows if it 
is noted that the Poincare section is defined using either g = or g = %. □ 

The lemma below is useful for finding fixed points of T p . 

Lemma 2.2. Let 

[•2irp 

<P P {lo,L ) = / Qidt, 
Jo 

where the arguments of Qi is evaluated with I = Iq + t/L\, g = go — t, L = Lq, and G = Gq. Then 
Mlo, (p/q) 1/3 ) = Mlo + 27T/p, (p/q) 1 / 3 ) and p (-/ o , (p/q) 1/3 ) = -<j> P {h, (p/q) 1/3 )- 

Proof. If we define 4>\(Iq,Lq) = f Q 2lT Vtidt, then L\ = Lq + (j)(lQ,Lo)^ + 0(fi 2 ) under the first 
return map T%. Using Lemma l2~Tl we may conclude that </>i(/o,^o) + <K — ^0 — 2tt/L 3 ) ,Lq) = 0. If 
Lq = (p/q) 1 / 3 , we have 

Mlo, (p/q) 1/3 ) + M-lo ~ 2nq/p, (p/q) l/3 ) = 0. (2.3) 

From T p = Tf, we get 

p-i 

(f) p (l ,L ) = J2Mlo + '2jn/L 3 ). (2.4) 
3=0 

That 4> p (Iq, (p/q) 1 / 3 ) has period equal to 2ir/p in Iq follows from Q2.4J1 and elementary number 
theory if it is noted that p and q are relatively prime. That 4> p (Iq, (p/q) 1 / 3 ) is an odd function of Iq 
follows from (f2~3|) . □ 
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3 Scaling and periodic points of the return map 



The variable A denned by L = (p/q) 1 / 3 + X^JJi can be used instead of L to blow up the region of 
the Poincare section near p/q resonance. The image of the point (Iq, Ao) in the l-X plane under the 
pth return map T p can be calculated using (|2.2|) . and it is given by 

h = lo ~ ciAo/x 1 / 2 + (c 2 X 2 + x (Z ))/i + 0(m 3/2 ) 

Ai = A + Hlo)v 1/2 + AoV'(^o)^ + 0(^ 3/2 ), (3.1) 

where 

d = faqWp- 1 / 3 and c 2 = 12vrg 5 /V 2/3 

/•2np 

<j>(l )= / Slidt 
Jo 

r2np r2np 3(7 4 / 3 f 2np 

x(io) = - n L dt-± n G dt-^— / siidtdr. (3.2) 

In (|3.2j) . the partial derivatives of Q must be evaluated at I = Iq + qt/p, g = go — t (where go 
is either or ir depending upon the choice of the Poincare section), L = (p/q) 1 / 3 , and G = 
(p/q) 1 / 3 ^ — e 2 ) 1//2 . The expression for T p given by 1)3.1 Jl and (|3.2|) is valid for any real value of lo 
and |Ao| < Sl/^ 1 / 2 - The domain of definition can therefore be taken as |Ao| < C\ with any positive 
constant C\ for sufficiently small \i. To derive the expression for ip(lo) given in 1)3.2)1 . we must 
use the second line of (|2.2j) and notice that the equation Lq = (p/q) 1 / 3 + Ao/u 1 / 2 and (|1.2|) imply 
Go = (p/q) 1 / 3 ^ — e 2 ) 1 / 2 + (q/p)XofJ 1 / 2 + 0(y). A term equal to 2-7rg has been dropped from the 
first line of ()3.1)1 as Iq and l\ are angles measured modulo 2ir. 

The function 0(Zo) equals 4> p (lo,Lo) defined by Lemma [2 . 1 1 when Lq = (p/q) 1 / 3 - By Lemma l2,ll 
<f>(lo) is an odd function with period equal to 2ir/p. Therefore, 4>(0) = and cp(ir/p) = 0. We make 
the following assumption about (ft(lo)'- 

Assumption A: For Iq £ [0,2-ir/p), the only points where 4>(lo) = are Iq = and 
lo = it /p. At these points, the derivative (f>'(lo) is nonzero. 

We turn to the verification of this assumption in Section 6. An example of a function which is 
odd with period 2n/p and which satisfies the assumption above is sm(plo). In fact, it will be shown 
later that 4>(lo) is proportional to sin(p/o)e' p_9 ' for small e. 

The assumption about 4>(Iq) can be put to use to find fixed points of T p in the l-X plane. We 
can use ()3.1)) and write 

(h - lo)/v l/2 = -ciX + (c 2 A 2 + x(lo))^ l/2 + O(ti) 
(Ai - A )/// /2 = </>(lo) + Xoi>(lo)n 1/2 + 0(ji). (3.3) 

When /x = 0, the right hand sides of the two equations in ()3.3)1 are both zero if Ao = and if lo is 
an integral multiple of ir/p. The implicit function theorem, along with the assumption about 4>(lo) 
stated above, allows us to infer that the right hand sides in 1)3.3)1 are for fi sufficiently small, 

= jir/p + 0(h), and A = (x(j^ Ip) / ci) fi 1/2 + Ofji), (3.4) 
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where j is an integer. Thus the points given by 1)3.4)1 are fixed points of T p for /j, sufficiently small. 
Since each application of the first return map T2 increments Iq by 2nq/p + 0(/u), where q and p 
are relatively prime, we may group the fixed points given by (|3.4|) into two sets, the first with 
j = 0, 2, . . . , 2(p — 1) and the second with j = 1,3,..., 2{p — 1) + 1. Then any fixed point in the 
first set moves to all other points in that set upon successive applications of Ti and returns to itself 
after the pih application; and likewise with the second set. The assumption about 4>(lo) implies 
that if one set of periodic points is elliptic then the other set is hyperbolic as will become clear 
shortly. 

If Assumption A about (f)(1) fails to hold, points Iq with (ft(lo) = and <P'(Iq) 7^ will still 
correspond to periodic points of 1)3.1)1 for /i sufficiently small. If Iq 7^ and £0 7^ 7Tj the corresponding 
periodic points can be the centers of asymmetric librations as described in the last section. 

4 Formal expansion of the stable manifold 

The expression for the map T p given by ()3.1|) and (|3.2j) can be rewritten by shifting the center of 
the l-X plane to (Z*,A*), where (Z*,A*) is the fixed point of T p given by 1)3 .4|) . with j being some 
integer. The map T p applied to the l-X plane centered at such a fixed point takes the form 

h = l - ciA /x 1/2 + (c 2 \l + x(h))n + r(l , A , ^)^ /2 

Ai = A + $(Z ),u 1/2 + A tf (Z )/i + *(*o> Ao, ^V 3/2 , (4.1) 

where 

X{h) = Xih+j^/p) - x(j'7r/p), $(i ) = <f>(l +jir/p), tf(Z ) = ^(k+jir/p). (4.2) 

The derivations of 1)2.2)1 and 1)3.1)1 imply that the remainder terms r and s in 1)4.1)1 are analytic in 
lo, Ao, and ^//I for any real Iq, Xq with |Ao| < C\, and sufficiently small in magnitude. Besides, 
r(0, 0, yfji) = s(0, 0, yfji) = because (0, 0) is a fixed point of 1)4.1)1 . 
The Jacobian <iT p of 1)4.1)1 at the origin is given by 

M*<(oU -t i/2 M x ' <0> * ( o))" +o c 3/2 )- ^ 3 > 

For // > and [i small, the fixed point is hyperbolic if $'(0) < and elliptic if $'(0) > 0. From 
(|4.2|) . it follows that $'(0) = (j)'(jir/p). The assumption of Section 3 implies that <^'(0) and 4>'(tt/p) 
are of opposite signs. Thus if the set of fixed points given by ()3.4|) is of hyperbolic or elliptic type 
for even j, the set of fixed points given by odd j must be of the opposite type. We shall assume that 
the fixed point used to shift the coordinate system and obtain 1)4.1)1 to be of hyperbolic type, which 
means $'(0) < 0. Let a = ^/-$'(0)/ci. Then a calculation using ()4.3)1 shows that the eigenvalues 
1 - c?an x l 2 + ((x'(O) + (0))/2)/x + 0(/i 3 / 2 ) and 1 + a 2 ^ 1 / 2 + (( x '(0) + *(0))/2)// + 0(^ 3 / 2 ) of 
cZTp correspond to eigenvectors of slopes 

a H fj, 1 +0( / u) and -a H \i 1 +0( / u), (4.4) 

ZCl ZCl 

respectively. The slope of the stable manifold of the fixed point of the map ()4.1)1 at the origin must 
be given by the first of the two expressions in 1)4.4)1 . When we derive an approximation to that 
stable manifold, 1)4.4)1 will serve to check the correctness of that approximation. 
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jc/p 



Figure 2: In the figure above, u(l) has been scaled down by 0.1 to make it fit. The plots correspond 
to the case q/p = 3/1, e = 0.1. 

To find the stable manifold of this fixed point, we rewrite (|4.1|) in the following form: 

Cl Cl 

\\ = A 2 + 2A <I>(/o)^ 1/2 + mo? + 2A§*(Z ))M + • • • (4-5) 

Let (Iq,Xq) be a point on the stable manifold. Its iterates (Zi, Ax), (h,^2), and so on are also on 
the stable manifold. In addition, l n — > and A n —* as n —* oo. Using (|4.5|) . we may write 

A 2 = A 2 , + 2(A $(/o) + Ai*(ii) + • • • + A n _i<I>(/ n _i))^ 1 / 2 + 0(/i) 

= A 2 , - -((/i - Z )$Go) + [h ~ hMh) + ■■■ + (In- *n-l)*(*n-l)) + 0(M). 

^1 



Note that 



Using (|4.6|) and noting that — lj is 0(/i 1//2 ), we have 

/fin 



O fin 

a 2 =a 2 --/ ^(Odz + o^). 



Taking the limit n — > oo, we find that formally the stable manifold is given by A = u{l) +0(/// 2 ), 
where 



u(l) 2 = -— f $(l)dl. 

Cl Jo 



(4.7) 



The positive root must be used if I > and the negative root if I < 0. It can be verified that this 
expression for u(l) agrees with (|4.4|) for the slope at the origin. Figure El plots u(l) and U(l) = u{l) 2 . 
To find the next term in the expansion of the stable manifold, we use ()4,5(l to get 

n— 1 n—1 

X 2 n = Xl + 2///2 £ A^) + M £ *(Z,) 2 + 2Aj*(i i ) + • • • (4.8) 
j=o i=o 

and use (|4.5|) and (|4.6j) to get 

Cl Cl 

= -- $(/) dZ + - (c 2 A 2 <I>'(^)/2 + (c 2 A 2 + x(ii))*(«i))M + • • • (4-9) 

^1 Cl 

Using (|4.8|) and (|4.9j) and by turning a sum into an integral as before, we get the expansion 
A = u(l) + v (l)n 1 / 2 + O(n) for the stable manifold , where 

= i /' 5« + gjCg+xW^o + pi + „(,)*(,) ( 4.io) 



ciu(Z) 7 2 v ' 2u(i) 

This formula for v(l) agrees with ()4.4|) with regard to the slope at the origin. This procedure can 
be repeated to calculate more terms in the expansion of the stable manifold. 



5 Homoclinic points near resonances 

Assume that the stable manifold of the origin under the map (|4.1|) is the graph of the function 
A = Mq(1). Then Mq{1) must satisfy a functional equation of the form 

A + $(ZoV 1/2 + • • • = M (l - Cl \ n 1/2 + •••)■ 

Assume that Mq(1) = u{l) + v{l)n 1 ^ 2 + w(l)n, where u(l) and v(l) are given by (|4.7|) and l|4.10|) . 
respectively. Then w(l) must satisfy the functional equation 

(1 + c lU '(lo)n 1/2 )w(lo) = w(h) + s 2 (/ , ^w(l ), ^n)n 1/2 , (5.1) 

where s 2 is analytic in its arguments for < Iq < 3ir/2p, ^JJiw real and bounded by a large constant, 
and L/jul < y^Tio for some no > 0. In addition, s 2 (0, 0, y/JL) = 0. In Lemma l^Tl below , we prove that 
(|5.1|) has a unique C 1 solution w(l), with < I < 37r/2p, for < \i < no and some no > 0. We also 
prove that \w(l)\ and |tu'(Z)| are bounded by constants which are independent of n but which may 
depend upon no- Therefore, the stable manifold of Q4.1J) is the graph of A = u(l) + v(l)n 1 ^ 2 +w(l)n 
for < / < 3n/2p and < /i < no, where u(l) and v(l) are defined by (|4.7|) and l|4.1U|) . 

As discussed in ^3], the functional equation (|5.1|) for w is obviously a contraction for < I < 
tt/p — e, e > 0. But the construction of homoclinic points requires the existence of w to be proved 
over a larger interval as in the lemma below. The proof of the lemma uses a technique found in 
|14j . Another approach can be found in 

Lemma 5.1. For some no > and any n £ (0, no}> there exists a unique C 1 function w(l) such that 
1)5. lj) is satisfied for < lo < 3ir/2p, with \ w(l)\ and \w'(l)\ bounded by constants for < I < 3ir/2p. 
The constants are independent of n but may depend upon no ■ 
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Proof. We will look for a continuous solution of (|5.1j) that satisfies w(Q) = 0, |exp(— Kl)w (Z)| < C*, 
and Lip(exp(— Kl)w{l)) < L* , where Lip(-) is the Lipschitz constant. The choice of the positive 
constants K, C*, and L* will be made later in the proof. 

The function r(Z, A, ^/Jl) from the first line of (|4.1|) and the function s 2 (^ v^ 10 ' V^) f rom t|5 . 1 f) 
determine the initial choice of hq. Let w* = fiw and consider r*(Z, w*, ^JJi) = r(l,u(l) + ^/JIv(l) + 
w* iy fji). The constant fiQ > is chosen so that r* (Z , w* , ^fjl) is analytic in its arguments over the 
compact domain D where < Z < 3ir/2p, w* is real and |io*| < C, and \^fji\ < \tJJiq\- By taking 
/io small enough, we can assume C to be as large as we please. It is enough to assume C to be 
twice the height of u(l) depicted in Figure |2 for example. Now let w* = ^[Jiw and assume the 
choice of /io to be such that s 2 (Z,w*, yfji) is also analytic in its arguments in the compact domain 
D. We note r* (0,0, ^//Z) = 8 2 (0,0,y/p) = °- 

The proof, which is organized into a number of steps, introduces many constants. The constants 
that depend on the domain D will be denoted by subscripting D. The constants that do not depend 
upon the domain D will be denoted by subscripting C. The constants that depend upon D are 
typically upper bounds for the magnitudes of derivatives of r* and s 2 over the domain D. The 
constant /to may be made smaller by some of the steps in the proof. But the bounds obtained using 
the domain D as specified above will of course apply even if /xo is made smaller. All constants 
introduced in the proof are strictly positive. 

1. For < Iq < 3ir/2p, by (|4.1j) Zi as a function of Iq is given by 

h = l Q - ciA /x 1/2 + (c 2 X 2 + X{h))v + r(l , A , Vt)v 3/2 , (5.2) 

where Ao = u(Iq) + u(Z )/j 1//2 + w(Iq)/j. It is possible to think of (|5.2|) as defining l\ in terms 
of Iq and w. By the assumption about <p(l) in Section 3 and (|4.7|) . it follows that u'(0) > 
and that u(l) > for < I < 2>Tr/2p. As u(0) = 0, there must be a constant C\ such that 
c\u{l) > C\l for < Z < 3ir/2p. Both the // and /z 3 / 2 terms in (|5.2|) vanish when Iq = and 
w = 0. Further, |u>(/o)| < Lip(w)?o- Thus the magnitudes of the two terms can be upper 
bounded by (C 2 + C3 Lip(u>)/i)Zo/i and (Z?i + D 2 Lip(u>) //)Zo/i 3//2 , respectively. Therefore, we 
may assert h < (l — (Ci/2)^ 1/,2 )io for < I < 3ir/2p and < /i < /io, with /iq made smaller 
if necessary. 

2. Let l\ be obtained using (|5.2|) with Zo replaced by Zq but with the same w. Assume < 
^Oi^o — 37r/2p and let C4 = Lip(citi(Z)). Then, as in the previous step, it follows that 
\h — l*\ < (l + 2C 4 /i 1//2 ) I Z — Zq I for < fi < fi and /x sufficiently small. 

An additional fact about \l\ — Z|| will be needed. Assume < Zq,Zq — K/2p. By the mean 
value theorem, ci(u(lo) — u(Iq)) = ciu'(1)(Io — Iq), and by the assumption about 4>(l) in Section 
3 and gSJ), u'{l) > C 5 for < I < n/2p. Therefore, \h - Z*| < (l - {C^/2)n l l 2 ) |/ - Z*| or 
simply \li — Z*| < \Iq — Zq | in this situation, for < fi < fiQ and /io sufficiently small. 

3. From (|5.1|) . we may obtain the following iteration: 

/ \ /, v P /, n / Tsi \ (1 \ 1 exp(-^o)g 2 (Zo,07w n (Z o ),^/7i) 1/2 
exp(-iiCZo)w n+ i(Zo) = F u , n (Z )exp(-KZi)u; n (Zi) H — ... r 1/2 /i 7 , 

(5.3) 

where l\ is obtained from Zo using 1)5.2(1 but with u; replaced by w n , where w n belongs to the 
class of functions for w specified at the beginning of this proof, and where the contraction 
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factor F Wn is given by 

exp(K(Zi - Z )) 

1 + CiU^/oJ/U 1 ^ 

The constant -fT will be chosen so as to make F Wn a sufficiently strong contraction. 

First consider < Iq < ir/2p. Let the minimum value of c\u'{Iq) for Iq in this range be Cq. 
Since h - l < -C 1 l [i l / 2 /2 by the first step, it follows that F Wn (l ) < (l - (C 6 /2)/i 1/2 ) for 
< Iq < vr/2p, < \i < fj,Q, and //o sufficiently small. 

Next consider 7r/2p < Zo < 37r/2p. Let the minimum value of c\u'{l) for Iq in this range be 
-C 7 . Then 

exp(-TO 7 r// 1 / 2 /4p) 
F Wn (lo)< 1 _ C7fil/2 • 

Choose K so that Kdn/Ap > 5C 4 + 2C 7 and conclude that F Wn (l ) < (1 - 4C 4 ^ 1 / 2 ) for 
7r/2p < /q < 37r/2p, < /j < /io, and /^o sufficiently small. 

4. By assumption, |exp(— K l)w n (J)\ < C* for < I < 3n/2p. Using (|5.3[) . we may upper bound 
|exp(-in)it; n+ i(0| by (1 - C 8 n 1/2 )C* + D^ 1 ' 2 , where C 8 = min(C 6 /2, 4C 4 ) and D 3 is an 
upper bound of the coefficient of /i 1//2 in the last term of (|5.3|) . By choosing C* > D 3 /Cs, we 
assert that |exp(— Kl)w n+ i{l)\ is also upper bounded by C*. 

5. By assumption, Lip(exp(— Kl)w n (l)) < L*. Let 

Q = ex.p(-Kl )w n+ i(l ) - exp(-KQw n+ i(lo), 
where < Zq < Iq < 3n/2p. We will upper bound \Q\. 

Using (|5.3j) . both the terms of Q can be replaced by expressions in terms of w n . The resulting 
expression for Q equals A\B\ — A2B2 — where A\ = F Wh (Iq), B\ = exp(— kl\)w n (l\), A2 = 
F Wti (Iq), and B2 = exp(—kli)w n (l*) — plus another term which equals the difference of two 
quantities times /x 1 / 2 . This other term will be denoted by Q r \j^l 2 . 

To bound \Q\, first consider the case Iq > ir/2p. We write |^4i-Bi — ^2-82! < \A\\ \B\ — £2 1 + 
I-B2I \A\ ~~ -^2 1- By the third step and the assumption about Iq, \A\\ < (l — 4C4/i 1//2 ), and 
\Bi — -B2I < L* \li — ll\ < L*(l + 2C4/U 1 / 2 ) \Iq — Iq\, where the last inequality follows from 
the second step above. Therefore, \Ai(B\ — B 2 )\ < (l — C4/j 1//2 )L* \l — Zq | . A simple esti- 
mate shows that ~Lip(F Wn (l)) < Cg^u 1 / 2 for sufficiently small hq. Therefore, | ^2(^-1 — ^2)! < 
C*C 9 M 1/2 \h -Iq\. To upper bound Q r , note that the coefficient of /j 1//2 in (|5,3j) has a Lipschitz 
constant with respect to Iq that can be bounded as D4 + L>5 Lip(w n )/j, 1 ^ 2 or D4 + D5L* /x 1 / 2 . 
Therefore 

101 < ((1 - CV/ 2 )L* + C^g/i 1 / 2 + D 4 /i 1/2 + A^>) Ko - • 

If L* > 2(C*C 9 + D 4 )/C 4 , then \Q\ < L* |Z - Zg| for < // < /j and /i sufficiently small. 

Consider the case < Iq < Iq < ir/2p. In this case, the argument is identical to that given in 
the previous paragraph, except that the bound on \Ai\ must be replaced by (l — (Ce/2)/i 1 / 2 ) 
form the third step, and \B\ — < L* \l\ — l*\ < L* \Iq — Iq\ from the additional fact in the 
second step. In this case, if L* > 4(C*C 9 + D 4 )/C 6 , then \Q\ < L* \l - Zq | for < fi < /jlq 
and /Iq sufficiently small. 

The choice L* = 2(C*C 9 + D 4 ) max(l/C 4 , 2/C 6 ) implies Up(exp(-Kl)w n+1 (l)) < L*. 
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6. If w n (0) = 0, then w n+ i(0) = since ^2(0, 0, yfp) = 0. This observation together with the 
choice of C* and L* in the fourth and fifth steps implies that w n +\ belongs to the same class of 
functions as w n . The third step with some other estimates given above implies that the map 
w n — ► w n+ i given by 1)5.3)1 is a contraction for < \i < /iQ. We conclude that there is a unique 
continuous solution w of 1)5.1)1 such that \exp(—Kl)w(l)\ < C* and Lip(exp(— Kl)w(l)) < L* 
for < I < 3vr/2p. 

7. If w(l) is continuously differentiable its derivative can be easily bounded in terms of K, C* , 
and L* . To complete the proof, it suffices to show that w(l) is continuously differentiable. 
The standard stable manifold theorem states that w(l) will be analytic in I in a neighborhood 
of I = 0. The stable manifold over the interval < / < 3ir/2p can be obtained by repeated 
applying T~ l to a local segment. Therefore w(l) must be continuously differentiable. 

□ 

If (j)(lo) defined by 1)3.2)) satisfies the assumption in Section 3, the fixed points of the map T p 
given by 1)3.4)1 are hyperbolic for j even or for j odd. One of these hyperbolic points was shifted 
to the origin in 1)4.1)1 . and we proved that the stable manifold of the origin is given by the graph of 
A = u(l) + v{l)^L 1 / 2 + 0(/i) for < I < 3ir/2p. In the L-l plane, the stable manifold is the graph of 

L = (?) V3 + ^ + U {1- jTr/p)^ 2 + v(l - jtt/pV + 0(/// 2 ) (5.4) 

\q/ ci 

for < I — jir/p + O(fi) < 3ir/2p, where \ is given by ()3.2|) . u is given by 1)4. 7|) . and u is given by 

The choice of the Poincare section as either g = or g = tt is yet to be made. To facilitate the 
construction of homoclinic points, it is also useful to pick j in ()3.4|) carefully. There are four cases. 

• If p is odd, then g = is chosen as the Poincare section. If (p'(0) > 0, then j = —1. 

• If p is odd and (j>'(0) < 0, then j = (p— 1). 

• If p is even, first try g = as the Poincare section. If (j)'(0) > 0, then j = —1. 

• If p is even and </>'(0) < with g = as the Poincare section, choose the Poincare section 
g = tt and j = —1. 

Theorem 5.2. Let p and q be relatively prime positive integers and let p/q ^ 1/1. Assume that 
e lies in the interval [e m i n , e max ] defined at the beginning of Section 2. Let <f)(lo) be defined by 
()3.2)) with L = (p/q) 1 / 3 and G = (p/q) l ^{^ — e 2 ). Assume that <p(lo) satisfies Assumption A of 
Section 3, namely, for l £ [Q,2n/p), 4>{Iq) = only if l = or l = n/p and 4>'(Iq) ^ at those 
two points. Identify the Poincare section for the flow of the Hamiltonian 1)1.2)) of the restricted 
three-body problem with a region of the L-l plane by using H = —(p/q)~ 2 / 3 /2 — {p/q) 1 / 3 ^ — e 2 ) 
and by choosing g = or g = tt as indicated above. Then the pth return map T p given by (|2.2[) has 
a homoclinic point on this Poincare section at (lh,Lh), where Lh = (p/q) 1 / 3 + u(ir /p)^ 1 / 2 + O(fi) 
and lh = or 1^ = tt, for < /i < /io and po sufficiently small. 

Proof. We give a proof for the second case listed above. The other cases are treated similarly. In 
this case, p is odd, the Poincare section is g = 0, and j = (p— 1). By Lemma 15. li the representation 
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of the stable manifold given by Q5.4j) is valid for < l — (p — l)ir/p+0((j,) < 3tt/2p, and therefore the 
stable manifold crosses the line I = tt. By Lemma 12.11 this stable manifold can be reflected about 
the line I = w to obtain an unstable manifold. Thus we find a homoclinic point with I = tt. □ 

The homoclinic point constructed in Theorem 15.21 can be mapped using the first return map T\ 
to obtain a ring of p homoclinic points with > (p/q) 1 / 3 . In Section 4, we constructed the stable 
manifold of the fixed point of (|4.1|) at the origin over < lo < 3ir /2p. A similar construction applies 
over the interval —2>ir/2p < Iq < 0. That construction can be used to find a ring of p homoclinic 
points with < (p/q) 1 / 3 - 

6 Verification of the condition on (f)(lo) 

Equations (|1,3|) and (|l,4j) define f2 as a function of L, I, G, g. If the angles / and g are replaced by 
—I and —g, 0, is unchanged. Therefore, O can be Fourier expanded as 

SI = ^c mn cos(m/ + ng), (6.1) 

m,n 

where m can be any non-negative integer and n can be any integer. The coefficients c mn are 
functions of L and G. We use L = (p/q) 1 / 3 and G = (p/^) 1//3 (l — e 2 ) 1 / 2 . The Fourier expansion 
is valid if 1 + r 2 — 2rcos(#) > and is therefore valid if e is sufficiently small for p/q ^ 1/1. The 
coefficients c mn can be expanded as power series in e, and it is possible to determine the precise 
radius of convergence of these series. For our purposes, it suffices to note that all these series 
converge in some neighborhood of e = 0. If O is differentiated with respect to L or G the condition 
e > has to be imposed. However, the partial derivatives of O with respect to / and g and O itself 
are analytic in a neighborhood of e = 0. 

Lemma 6.1. // the Fourier coefficients c mn of (|6.1|) are expanded in powers of e, the lowest power 
of e with a possibly nonzero coefficient is e\ m ~ n \ . 

Proof. The quantity £1 can be expanded as Ylm=o c m( r ) cos(m#). First consider e = 0. Then 
r = (p/q) 1 / 3 and 9 = v + g = l + g. Therefore the only nonzero terms in the expansion (|6.1|) occur 
when m = n. If e ^ 0, then r depends upon I and v is no longer equal to I. This lemma follows 
when the dependence of r and v on I is taken into account. The way to do this can be found on 
pages 44 and 170 of ^01 or on page 35 of □ 

By (|3.2|) . 4>(Iq) = Jg 27rp &i dt, where Qi must be evaluated at / = lo+qt/p, g = go—t, L = (p/q) 1 / 3 , 
and G = (p/q) 1 / 3 ^ — e 2 ) 1 / 2 ; go = or go = it depending upon the choice of the Poincare section. 
Using (|6.1|) . we get 

oo 

(f)(l ) = -2-Kp 2 mc mp<mg sm(kpl + kqg ). 

m=l 

Let c p>q = c*(p,q)e\ p - q \+0(e\ p - q \ +1 ). Then hv LemmalO cj)(l ) = ±2irp 2 c*(p, q) sin(p/ )e |p - g| + - • • . 
Thus the assumption about <p(lo) in Section 3 and in Theorem 15.21 will be verified for e > and 
e small if we can show that c*(p,q) ^ 0. An expression for c*(p,q) can be obtained from the 
discussion of the quantity 

67rg 4/3 f 2np 

C(e,p,q) = T-pr- I ^udt 

P 1 ' 3 Jo 
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Figure 3: Plots of (f)(1) against I with q/p = 3/1 and with e ranging from 0.1 to 0.8. 
given in If p < q, for example, 



evaluated at a = (p/q) 1 ^ 3 ; above D stands for the differential operator a-^ and b n (a) are defined 
by the expansion 



The b n (a) are hypergeometric functions whose series converge for |a| < 1. The value of c*(p,q) can 
be obtained using the expression given above or by other means. 

In Figure E3 we have plotted <p(l) with q/p = 3/1. From that figure, it is clear that the 
assumption about 4>(Iq) is valid for even large values of e. 

7 Resonance boundaries 

Discussion of the averaged circular restricted problem and its use in sketching the boundaries of 
q/p Kirkwood gaps in the a-e plane can be found in [Hj, 0, |Hj? El and [Ej- The basic procedure 
is to average the Hamiltonian (J1.2|) by retaining only the terms in the Fourier expansion (|6.1|) of 
Q with ml + ng = k(pl — qg) for some integer k. These are the resonant terms. In some instances 
such as certain additional terms are added to model the effect of the eccentricity of Jupiter 
and the secular variation of its elements. 

This averaged Hamiltonian has 1 degree of freedom. Its fixed points and separatrices are used 
to approximate the boundaries of resonance. For small values of the asteroid eccentricity e, some 
of the nonresonant terms dropped during averaging have larger coefficients than any of the terms 
retained during averaging. To some extent the influence of the nonresonant terms is captured by 



c*(p, q) 




(ab q (a)) 




n=— oo 
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q/p 


3/2 


2/1 


7/3 


5/2 


3/1 


4/1 


e 


.10 


.13 


.08 


.08 


.07 


.09 



Table 1: This table gives the values of e for certain q/p below which the periodic points of Section 
3 fail to exist for p, = 10~ 3 . Boundaries of resonance obtained using the averaged circular restricted 
problem will not be valid below these values. 



q/p 


1/7 


1/6 


1/5 


1/4 


1/3 


1/2 


e 


.365900 


.320133 


.265532 


.199749 


.121094 


.036083 



Table 2: This table gives the values of e above which asymmetric librations exist in the unaveraged 
circular restricted problem for certain exterior resonances. 

the formal change of variables used to average the Hamiltonian, but this change of variables is often 
not taken into account. Even if it is, the averaging will not be valid at small values of e. 

The return map (|3.1|) corresponds to the unaveraged circular restricted problem. Only the 
resonant terms of £1 contribute to </>(/) defined by ()3.2|) . but the nonresonant terms contribute to 
both ip(l) and x(0- As e approaches 0, the magnitude of (f)(1) becomes much smaller than that of 
the other two functions in (|3.2|) . Thus for fixed /x and small e the periodic points (|3.4j) will not 
exist. 

These periodic points, when they exist, correspond to the fixed points of the averaged Hamilto- 
nian. The separatrices of the averaged Hamiltonian corresponds to the stable manifolds discussed 
in Sections 4 and 5. If for certain values of n,e,p,q the unaveraged circular restricted problem 
does not have the periodic points given by (|3.4|) . the fixed points and separatrices of the averaged 
Hamiltonian must be treated as artifacts of the averaging procedure. 

In Tabled we have given the minimum values of e required for the periodic points ()3.4|) to exist 
for some of the commonly studied resonances in the asteroid belt. The value of /u used is close to 
that of Jupiter. The boundaries of resonance obtained by averaging can be valid only above these 
values of the asteroid eccentricity. 

8 Asymmetric librations 

If the assumption about (j)(l) in Section 3 holds, the circular restricted problem has two resonant 
periodic solutions for [i small. One of these is of elliptic type and therefore there will be solutions 
that librate around the periodic points that correspond to it in the Poincare section given by the 
l-L plane. These are symmetric librations. 

As shown in Figure 01 the assumption about (f)(1) in Section 3 can fail for some exterior reso- 
nances as e increases. In both the plots shown in that figure, the periodic point with / ~ ir/p is 
initially of elliptic type as the slope of (f)(1) is positive, and when e increases, it undergoes a pitchfork 
bifurcation and becomes a point of hyperbolic type. Elliptic points appear on the Poincare section 
at values of I that are not 0([i) close to any integral multiple of vr/p. The librations around these 
points are termed asymmetric. 

A study of asymmetric librations in the exterior resonances using averaged equations can be 
found in [2]. In Table |H1 we have listed values of e above which asymmetric librations occur in the 
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Figure 4: The two figures show plots of (f)(1) vs. I for I € [0, 2ir/p) for q/p = 1/3 and q/p = 1/7, 
respectively. The values of e are 0.08, 0.12 and 0.16 in the left plot, and .30, .36 and .40 in the 
right plot. 

unaveraged equations for some exterior resonances. If the averaged equations imply the existence 
of asymmetric librations below these values of e, those must be considered artifacts of averaging. 
However, if e exceeds these values, asymmetric librations will be found for small enough. 
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